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红 精 灵 发 光 光 谱 的 数值 模拟 研究 


江 F X857 KER’? ZR? Ium 
1( 中 国 科学 院 空间 科学 与 应 用 研究 中 心 ”北京 100190) 
2( 南 京 航空 航天 大 学 航空 宇航 学 院 ” 南京 210016) 


摘 要 红 精 灵 是 发 生 在 雷暴 层 云顶 的 一 类 大 气 瞬 态 发 光 现 象 ， 是 能 量 由 对 流 层 耦合 到 中 高 层 大 气 的 直接 证 据 . 
其 发 光 光 谱 研 究 是 了 解 整个 事件 对 中 高 层 大 气 能 量 注 入 的 重要 手段 ， 有 助 于 认识 事件 发 生 区 域 的 大 气 电离 度 及 事 
件 过 程 的 能 量 电子 分 布 , 进一步 为 研究 红 精 灵 的 产生 机 制 提供 重要 信息 , 同时 为 大 气 辐射 背景 资料 研究 提供 重要 依 
dE. 本 文 利用 Boltzmann 方程 求解 了 电场 作用 下 弱电 离 气体 中 电子 能 量 分 布 的 时 变 函 数 ， 以 此 为 基础 , 模拟 计算 
了 红 精 灵 各 典型 发 射 带 的 光谱 强度 . 模拟 计算 结果 表明 , 约 化 电场 E/N 越 强 , 电子 获得 的 能 量 就 越 多 , 高 能 电子 也 
就 越 多 , 致使 撞击 中 性 大 气 产生 的 辐射 光 强 就 越 强 ; 模拟 显示 红 精 灵光 辐射 谱 分 布 从 远 紫 外 直至 近 红 外 . 

关键 词 红 精 灵 , 大 气 瞬 态 发 光 现 象 , 能 量 分 布 , 体 发 射 率 
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Abstract Sprites are the transient luminous events above very active thunderstorms and are the 
direct evidence of the energy in troposphere coupled to the upper atmosphere. The emissive spec- 
trum study is the important means to understand the energy injected into upper atmosphere and 
regional ionization, the electron energy distribution during the whole event. These parameters are 
useful for the study of sprites’ physical mechanism, in addition as a kind of natural atmospheric 
luminous phenomenon, the emissive spectrum study is useful for providing important information to 
the research of atmospheric radiation background information. Based on the Boltzmann Equation 
to calculate the time evolution of the electron energy distribution function in a mixture of partially 


ionized gases with an applied electric field, the paper calculated the spectral intensity in typical 
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emissive band for sprites. Simulation results show that the emissive intensity is more intensive with 


more intensive E/N and the emissive spectrum is from far ultraviolet to near infrared. 


Key words Sprite, Transient luminous events, Energy distribution, Volume emissive rates 


0 引言 


红 精 灵 是 1989 年 Franz 等 Ul 在 一 次 低 光 度 相 
机 测试 中 意外 记录 到 的 瞬 态 大 气 发 光 现象 . 1994 年 
Sentman 等 外 利用 飞机 上 的 多 波段 相机 对 红 精 灵 进 
行 观测 , 发 现 红 精灵 发 生 的 高 度 在 40~90km, 上 下 范 
围 约 数 十 公里 , 通常 在 云 对 地 闪电 后 几 毫 秒 后 发 生 ， 
发 光 持 续 时 间 从 几 训 秒 至 数 百 毫秒 不 等 . 

红色 精灵 由 于 亮度 较 高 而 容易 被 观测 到 ， 成 为 
中 高 层 大 气 瞬 态 发 光 事 件 中 研究 最 多 的 现象 ， 准 静 
电场 模型 对 其 发 生机 制 的 解释 认为 ,闪电 放电 
后 雷暴 云 中 剩余 电荷 及 大 气 中 的 诱导 电荷 与 自由 电 
荷 在 中 高 层 大 气 中 产生 一 个 很 强 的 准 静 电场 (Quasi- 
Electrostatic Field, QEF), 并 渗透 到 中 间 层 和 低 电 高 
层 , 强度 随 高 度 和 时 间 变 化 并 产生 电子 加 速 , 进而 引 
起 大 气 的 击 穿 和 光 辐 射 . 1996 Æ Pasko 等 91 利用 经 
验 公 式 计 算 了 低 电离 层 中 的 电子 在 QEF 作用 下 引起 
的 电离 和 分 子 的 光学 激发 , 但 计算 是 通过 实验 数据 拟 
合 而 成 的 , 并 没有 考虑 电子 在 QEF 中 的 动力 学 过 程 . 
在 此 基础 上 , 本 文 考虑 电子 在 QEF 中 的 动力 学 行为 ， 
利用 Boltzmann 方程 计算 了 电场 作用 下 电子 与 中 性 
大 气 成 分 碰撞 的 动力 学 方程 , 进而 获得 电子 能 量 分 布 
函数 , 并 在 此 基础 上 计算 了 No 各 典型 发 射 带 系 的 发 
射 光谱 强度 . 


1 动力 学 模型 


部 分 电离 气体 中 存在 电场 时 ， 
速 作 用 使 能 量 增 大 , 与 气体 分 子 发 生 碰 撞 后 , 将 能 
传递 给 气体 分 子 , 造成 能 量 损失 ， R 态 
等 离子 体 放电 , 其 电子 能 量 分 布 人 化 不 再 服 
从 Maxwell 分 布 和 Druyvesteyn 分 布 Ul, 或 者 说 不 存 
在 普遍 适用 的 分 布 函 数 . 因此 , 需要 根据 Boltzman- 
n 输 运 方 程 计算 出 准确 的 电子 能 量 分 布 函数 , 此 方程 
可 以 描述 受到 电场 加 速 的 电子 与 气体 分 子 间 的 相互 
EH, 主要 考虑 二 体 碰撞 , 适用 于 电离 度 极 低 的 低 电 
离 层 . 


Boltzmann 输 运 方程 表示 如 下 : 


($ AME Zy DIGDD = (5). 


E = Eo sinwt. (2) 


其 中 , r 为 位 置 参数 , v 为 电子 速度 , e PT utar, 
m 为 电子 质量 , E 为 外 加 电场 强度 , w 为 电场 激发 频 
率 , Vr 和 Vo 分 别 为 空间 和 速度 梯度 算 符 . XX (1) 
中 f(7,v,t) 描述 的 是 电子 数 密度 在 空间 和 速度 上 的 
分 布 概率 , 其 函数 即 为 要 求 的 解 , 式 (1) 等 号 右边 是 
碰撞 过 程 中 影响 f 的 碰撞 项 . 

假设 电子 数 密度 在 空间 上 均匀 分 布 , 即 与 位 置 7 
无 关 , 那么 fr, vs t) > f(o 0. 将 (ost) 做 两 项 球 谐 


展开 (two-term spherical harmonic expansion), 有 
v 
f(t,v) = folt, v) + zf v). (3) 


其 中 , fo(v. t) 表示 球 对 称 的 电子 速度 分 布 函数 , 与 速 
度 大 小 有 关 ; fito) 表示 非 对 称 项 , EBAN 十 z 
轴 方 向 运动 , 此 项 可 写 为 filt, v) cos0, 这 里 9 为 速 
E v 与 电场 E RKA. 
考虑 气体 与 电子 间 动 量 的 转移 、 非 弹性 、 弹 性 
碰撞 及 超 弹 性 磁 撞 (电子 与 激发 态 气体 的 碰撞 ) 过 程 ， 
Boltzmann 方程 可 以 表达 为 如 下 电子 数 密度 的 形式 : 
On Ot Ja | 


Qt c Oe 


2 (et eane spe) 
— >_Ni. 
[Ri (E)n(e) - Ryle — esp)n(e = es. (9) 
其 中 ， 
2N2(E/N)?e n n), 
dn m(2e/m) MEE Leto Oa) (5) 


1/2 2 
Jas -n(2) > ös T- oele) 


ino 芳 等 : 红 精 灵 发 光 光 谱 的 数值 模拟 研究 


P-e) -erg (6 
Ryle) = (Z) ogl); U 
RO = (4) [Eten O 
N = EN? + NÌ); (9) 


n(&) = ne fole)? 
使 
| n(e)de = ne; 
(10) 
十 co 
E= | Jue ye de. (11) 
0 


X (1)~(11) rf, s 为 受 激发 气体 种 类 , 7 为 受 激发 气 
体能 态 , jf 为 从 电场 中 获得 的 能 量 , jo 为 弹性 碰撞 损 
失 的 能 量 , Pus; 为 电子 激发 s 气体 从 基态 到 激发 态 J 
的 碰撞 截面 , Ri 为 电子 与 激发 态 j 的 s 气体 产生 超 
弹性 碰撞 的 碰撞 截面 ,= 为 电子 平均 能 量 ，V 为 空气 

式 (4) 右边 第 一 项 为 电子 能 谱 的 基本 展开 , 为 简 
化 可 认为 放电 前 后 空气 数 密度 近似 保持 不 变 , 因此 通 
过 式 (4) 求解 可 获得 时 变 分 布 函 数 . 将 外 加 电场 转化 
为 B/N 的 形式 作为 输入 参数 , 这 样 处 理 的 目的 是 , 同 
时 兼顾 电场 的 作用 和 中 性 成 分 的 影响 , 且 减 少 了 方程 
的 输入 参数 . 以 E/N 为 输入 量 , 通过 有 限 差分 法 对 
方程 离散 化 , 反复 迭代 直到 电子 能 量 分 布 函数 达到 稳 
态 为 止 , 此 即 电子 能 量 分 布 函数 , 根据 式 (11) 可 计算 
电子 平均 能 量 . 

计算 中 假设 空气 是 N. 和 Oo 按照 体积 分 数 
为 2196 和 7996 的 比例 混合 而 成 , 除 N2 和 O2 的 
反应 外 , 其 他 粒子 的 碰撞 反映 情况 暂 忽 略 . 主要 考虑 
的 碰撞 过 程 及 其 对 应 的 能 量 阔 值 列 于 表 1 和 表 2, 各 
种 碰撞 截面 取 自 实验 结果 文献 [57-9]. 

图 1 给 出 的 是 当 E/N = 10Td (1Td = 
10-1V:cm2) 时 , 迭代 过 程 中 (达到 稳定 前 ) 电子 能 
量 分 布 函数 随时 间 的 变化 情况 .从 图 1 可 以 看 出 ， 
有 外 部 电场 存在 的 非 平 衡 态 下 ， 电 子 能 量 分 布 呈 明 
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显 的 非 Maxwell 分 布 ; 随 着 时 间 的 增加 , 高 能 区 的 电 
子 越 来 越 少 , 直至 达到 稳定 后 , 电子 能 量 分 布 函数 趋 
于 Maxwell 分 布 . 

图 2 为 1.0ns 时 刻 ( 非 稳定 时 刻 ) 电子 能 量 分 布 
函数 随 约 化 电场 (B/N = 10, 30, 50, 100, 200, 400, 
600 Td) 的 变化 情况 .从 图 2 可 以 看 出 , 在 E/N = 
10Td 时 , 大 部 分 电子 的 能 量 较 低 , 高 能 电子 比例 很 
少 , 电子 能 量 高 于 10eV 的 分 布 概率 几乎 为 零 ; 随 着 
约 化 电场 不 断 增强 , 最 大 分 布 概率 降低 , 能 量 分 布 范 
围 显著 增 大 , 低能 电子 减少 , 高 能 电子 增多 , 分 布 概 
率 向 高 能 方向 移动 , 越 来 越 背 离 Maxwell 分 布 . 

与 图 2 相对 应 , 图 3 给 出 了 同一 时 刻 B/N 与 电 
子平 均 能 量 的 关系 , 随 着 B/N 的 增 大 , 高 能 量 电子 
越 来 越 多 , 因而 平均 能 量 随 B/N 的 增 大 明显 增 大 . 


10° ey 
g ` 
$ 10? rr 
o ya 
8 x deest o a nme 
8 10 4 s 5 
E 106 
5 
$ 10° — t=0ns 
Pa zi 
5 10  t-0.05ns 
id 19 --t-0.5ns 
Ped 1072 — t23.0ns : 
[5] 
5 N2 E/N-10 Td, P-1Pa 


energy/eV 


图 1 电子 能 量 分 布 函数 随时 间 的 变化 (达到 稳定 前 ) 


Fig.1 Time evolution of the electron energy 


distribution function (before reach steady state) 


energy distribution function 


energy/eV 


图 2 ”电子 能 量 分 布 函数 随 约 化 电场 的 变化 (达到 稳定 前 ) 
Fig.2 Evolution of the electron energy distribution 


function with E/N (before reach steady state) 
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表 1 电子 与 氮气 分 子 的 主要 碰撞 过 程 


Table 1 Inelastic collisions of electrons with molecular nitrogen 
collision process reaction threshold energy /eV 
N2 rotational e 十 Na — e+N2 (rotational) 0.02 
Nə vibrational e4-N» Na2( —1 0.29 

e4N5 一 e+N2(v = 1 0.291 

e+N2 一 e4-N»(v = 2 0.590 

e4N» 一 e4-N»(v = 3 0.88 

e+N2 一 e-N»(v = 4 1.17 

e+N2 一 e+N2(v = 5 1.47 

e+N2 一 e+N2(v = 6 1.76 

e+N2 > e+N2(v = 7 2.06 

e+N2 N2(v=8 2.35 

Nə electronic e+N2 > e-No( A3 ET, v = 1-4) 6.17 
e+N2 — e+N2(A? E, v = 5~9) 7.00 

e+N2 — e-N3(B? II) [LPN] 7.35 

e+N2 Nə (W° A,) 7.36 

e+N2 一 e+N>2 (A't, v z 10) 7.80 

e+N2 一 e-N3(B? X7 ) 8.16 

e+N2 一 e+N2(a" X; ) 8.40 

e+N2 — e+N2(a' IT;) [LBH N;] 8.55 

e+N2 一 e+N2(w! Au) 8.89 

e+N2 — e+N2(C° Mu) [2PN?] 11.03 

e+N2 > e+N2(E’° XF) 11.88 

e+N2 > e+N2(a' Ef) 11.25 

N2 sum of singlet states e+N2 一 et+N+ 13.00 
Na ionization e+N2 2 2e-Ni (X? X7 , A? II) 15.6 
e+N2 — 2e+ Ni (B? XZ) [INNi] 18.8 


2 激发 率 的 计算 


得 到 电子 能 基 分 布 和 碰撞 截面 , 即 可 计算 得 到 各 
个 碰撞 过 程 的 激发 率 , 有 


十 co 
Kez k = Vk/N = n v2e/mef(e&)ex(e)de. (11) 


其 中 , vx (单位 s 1) 为 碰撞 频率 , 反映 的 是 一 个 电子 
在 有 电场 情况 下 碰撞 大 气 分 子 使 之 处 于 k 激发 态 的 
REJI; s 和 me 分 别 为 电子 能 量 (单位 eV) 和 质量 ( 单 


位 g), N 为 中 性 分 子 数 密度 ; f(e) 为 归 一 化 电子 能 量 
分 布 , 满足 hi 
MESE 


0 
ex (€) 为 各 反应 过 程 对 应 的 碰撞 截面 . 观测 和 数据 分 
WRH, 中 高 层 闪 电 辐 射 主要 源 自 No 或 Nz 的 激发 
辐射 , 而 Oo 分 子 的 激发 率 较 No 或 N# 低 得 多 P, 因 
此 本 文 未 予以 考虑 , 这 里 主要 考虑 1IPNa, 2PN;, LBH 
N2, Meinel N7, INN 5 个 带 的 激发 辐射 , 对 应 的 分 
子 或 离子 的 能 级 路 迁 及 激发 能 量 阔 值 列 于 表 300). 
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表 2 ”电子 与 氧气 分 子 的 主要 碰撞 过 程 


Table 2 


collision process 


O» rotational 


Inelastic collison of electrons with molecular oxygen 


e+O2 — e-4-O» (rotational) 


O» vibrational e4 0» 


e+02 一 e+02(v = 


e+02 一 e+02(v = 


O» electronic e+02 一 e--O»5 at Ag) 


O» dissociation 


eT O0» — e+ 


e+02 > e4 


e+02 — e-O(? P) + O(P) 


e+02 — e+O(3P) + O(? S?) 


O» ionization 


O» two-body attachment 


e+O2 — 2e+O# (X? II) 


reaction threshold energy /eV 
0.02 
> e+02 (v = 1) 0.19 
2) 0.38 
e+02 一 e+02(v = 3) 0.57 
4) 0.75 
0.977 
e+O2 一 e+O2(b1 57) 1.627 
e+02 — e4O»(c! X; ) 4.50 
6.00 
OGPP) + O(! D) 8.40 
-O(! D) + O(1D) 10.00 
14.7 
12.06 


e+ O2 一 e+05 = 


表 3 Sprite 各 谱 段 的 跃迁 能 级 、 激 发 阅 值 、 波 长 范围 及 与 中 性 大 气 碰撞 的 济 灭 率 


Table 3 List of major emission band systems and key parameters of molecular nitrogen 


WERE ”波长 /nm 


kana / (cm? 7) — kao, / (cms 7) 


N(B? II,) > Na (A? X) 


发 射 带 跃迁 能 级 
1PN2 7.35 
2PN2 N2 (C? Iu) > N2(B° Ig) 11 
LBH N2 Na(a! Ig) > Na(X? XF) 8.6 
Meinel N} N (A’ Xt) —> NI(X?5,) 16.7 
INN N (BP XT) > Ni (X° X7) 19 


图 4 给 出 了 激发 态 No(B?IL), Na(C° Mu), 
Na(a!I), NF (Æ$), Ni(BiEf) 对 应 的 激发 
Z kou 随 E/N 的 变化 . 由 图 4 UAH, kerk 
随 E/N 的 增 大 而 增 大 , 如 表 3 所 示 , 相对 前 三 个 No 
分 子 带 的 激发 ， 两 个 NI 离子 带 的 激发 阐 值 要 高 ， 
因此 当 E/N 较 小 时 ， 几 乎 没有 离子 带 的 激发 ， 随 
着 E/N 的 增 大 , 两 个 离子 带 的 激发 率 迅 速 增 大 , É 
至 超过 分 子 带 的 激发 率 . 


478~2531 1.6 x 1071 1.5 x 10-1? 
268-2546 1.42 x 1071 2.85 x 10-19 
100-260 23x10-H 4.3 x 10-1? 
550-1770 1.12 x 1071 2.9 x 10-1? 
286-.587 4.53 x 10-1? 7.36 x 10719 


3 体 发 射 率 计算 


高 能 级 态 的 分 子 / 离 子 向 低能 级 态 跃 迁 时 会 释放 
光子 , 光子 对 应 的 波长 和 由 两 个 能 级 态 的 能 量 差 决 
定 , 即 

A=1/(E,— Ei). (12) 
其 中 , ,为 高 能 级 态 的 能 量 , Ex 为 低能 级 态 的 能 量 ， 
单位 均 为 波 数 cm 71. 
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electron mean energy /eV 
[ej 


4 

2 

0 100 200 300 400 500 600 
(E/N)/Td 


图 3 E/N 与 电子 平均 能 量 的 关系 (达到 稳定 前 ) 


Fig.3 Time evolution of electron mean energy (before 


reach steady state) 


0 100 200 300 400 500 600 
(E/N)/ Td 
图 4 各 激发 态 对 应 的 激发 率 kerr B8 E/N 的 变化 
Fig.4  Excitation rate for the upper state as 
a function of E/N 


一 个 双 原 子 分 子 /离子 的 能 量 由 所 在 的 电子 能 
级 、 振 动能 级 及 转动 能 级 三 部 分 能 量 构 成 , 处 于 激发 
态 的 分 子 / 离 子 在 向 低能 级 跃迁 过 程 中 , 除了 电子 能 
级 的 变化 , 还 会 产生 振动 能 级 以 及 转动 能 级 的 变化 ， 
因此 不 同 的 能 量变 化 会 形成 不 同 波长 的 光 辐 射 ， 已 
知 某 激发 态 的 激发 率 korr, 可 以 计算 处 于 激发 态 的 
氮 分 子 / 离 子 数 密度 , 有 


Dk,» = NN, NeKez,k/ (TE, + kana NN + kgos Nos). 

(13) 
HH, Nw No, Ne 分 别 为 大 气 中 和 氮气、 氧气 及 
电子 的 数 密度 ; ka, OM k 激发 态 氮 分 子 / 离 子 与 大 
气 中 Na 碰撞 引起 的 滩 灭 率 (cm3.s-!)，kso。 为 激 
发 态 氮 分 子 /离子 与 大 气 中 氧 分 子 Os. 碰撞 引起 的 
淳 灭 率 (cm3.s-!)， 具 体 数值 如 表 1 所 列 ; mw = 


Chin. J. Space Sci. 空间 科学 学 报 2015, 35(3) 


(SAS) 为 大 激发 态 、 v 振动 态 的 辐射 寿命 Gn 


kv 
位 s), 这 里 ATE 为 从 高 能 态 到 低能 态 跃 迁 的 Einstein 
系数 (单位 s D). 
已 知 激发 态 的 氮 分 子 /离子 数 密度 nkv， 可 以 计 
FH k MES, v 振动 态 跃 迁 到 Kr WES, v 振动 
态 的 体 发 射 率 , 即 


DEK (A) = np qt ABE, (14) 


其 中 , aor 为 从 基态 x 的 第 0 振动 态 激 发 到 激发 态 
的 第 v 振动 态 的 Franck-Condon 系数 H. 
将 式 (13) 代入 式 (14), 可 以 得 到 


kk' xk Akk' 
e (4) =Nk, vlov Ay, = 
xk Akk' 
Ny, Nekez kqov Ay, / 


(ris 十 kw; Ny, + kq02 No,). 
(15) 


对 于 中 高 层 瞬 态 发 光 事 件 发 生 的 高 度 区 域 ， 大 
气 氮气 、 氧 气 和 电子 的 数 密度 基本 不 受 太阳 活动 以 
及 时 空 变化 的 影响 ， 那 么 对 于 某 个 发 射 带 而 言 ， 体 
发 射 率 则 就 主要 取决 于 激发 率 kex 图 5~7 计算 
了 80km 高 度 处 E/N = 20, 100, 200, 400 Td Ff, 
各 发 射 带 的 体 发 射 率 随 波 长 的 分 布 . 由 图 5~7 可 以 
看 出 , 各 发 射 带 的 体 发 射 率 随 E/N 的 增 大 而 增 大 ， 
在 E/N 较 低 如 小 于 100 Tq 时, 光 辐 射 主要 以 几 个 
分 子 发 射 带 为 主 , ME E/N 的 增 大 , 两 个 离子 发 射 
带 的 体 发 射 率 逐 渐 增 强 , 量 级 与 分 子 发 射 带 的 体 发 身 
的 光谱 分 布 从 远 紫 外 至 近 红 外 . 


4 结论 

根据 Boltzmann 动力 学 方程 的 数值 求解 计算 出 
以 E/N 为 函数 的 电子 能 量 分 布 , 通过 电子 能 量 分 布 
与 大 气 成 分 的 激发 截面 进一步 得 到 以 B/N 为 函数 的 
激发 率 系数 , 在 此 基础 上 , 模拟 计算 了 红 精 灵 各 典型 
发 射 带 的 光谱 强度 . 结果 表明 , 约 化 电场 B/N 与 发 
射 带 的 光 强 存在 定量 关系 , 约 化 电场 E/N 越 强 , 电 
子 获得 的 能 量 则 越 多 , 从 而 导致 高 能 电子 越 多 , 撞击 
中 性 成 分 产生 的 光 辐 射 就 越 强 . 模拟 的 结果 表明 , 红 
精灵 光谱 范围 从 远 紫外 直至 近 红外 . 
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图 5 发 射 带 1PN2 和 2PN» 在 B/N = 20, 100, 200, 400 Td 时 各 发 射 带 的 体 发 射 率 随 波 长 的 分 布 
Fig.5 Volume emission of emissive band 1PN2 and 2PN» at E/N = 20, 100, 200, 400 Td 
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图 6 发 射 带 LBHN2 和 Meinel N7 在 E/N = 20, 100, 200, 400 Td 时 各 发 射 带 的 体 发 射 率 随 波 长 的 分 布 
Fig.6 Volume emission of emissive band LBH N» and Meinel Ni at E/N = 20, 100, 200, 400 Td 
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图 7 发 射 带 INN3 在 E/N = 20, 100, 200, 400 Td 时 各 发 射 带 的 体 发 射 率 随 波 长 的 分 布 
Fig.7 Volume emission of emissive band 1NN+ at E/N = 20, 100, 200, 400 Td 
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